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We investigate the compositional dependence of the total energy of the mixed 
crystals (Ga,Mn)As co-doped with As, Sn, and Zn. Using the ah initio LMTO-CPA 
method we find a correlation between the incorporation of acceptors (Mn, Zn) and 
donors (Sn, antisite As). In particular, the formation energy of AsQa is reduced 
by approx. 0.1 eV in the presence of Mn, and vice versa. This leads to the self- 
compensating behavior of (Ga,Mn)As. 

PACS numbers: 71.15.Ap, 71.20.Nr, 71.55.Eq, 75.50.Pp 

I. INTRODUCTION 

The Mn doping of the III-V semiconductors has two principal effects in the III-V diluted 
magnetic semiconductors (DMS). The correlated d-electrons at Mn atoms form local mag- 
netic moments and the hybridization of the d-states with the band states results in various 

n 

magnetoelectric and magnetooptical phenomena [1| . In addition, Mn atoms substituted 
for a trivalent cation act as acceptors and introduce holes into the valence band. It is now 
generally accepted [J] that the ferromagnetic coupling between the local moments in the 
III-V DMS is mediated by the mobile holes in the valence band. 

In reality, however, the number of the holes is much smaller than the nominal concentra- 
tion of Mn 0,0] This indicates that a large amount of compensating donors is present, with 
a strongeffect on both conductivity and Curie temperature of these materials. As pointed 
out in p, the nearly constant (or even decreasing) doping efficiency of order 0.1 ~ 0.2 
can be explained only assuming that the number of the donors increases proportionally to 
the concentration of Mn. This was the reason to suggest that Mn atoms in the interstitial 
positions - being double donors - might have a principal role in the compensation Q . 
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Another favorite candidate for the compensating donor is an As antisite defect. These 
defects are well known in crystals grown with an excess As, in particular in the (Ga,Mn)As 
films 0]. So far, however, there was no reason to expect any correlation between the con- 
centration of As antisites and the concentration of Mn. 

Only recently, an increase of the number of the As antisites was indicated by comparing 
Curie temperature obtained from ab initio calculations for (Ga,Mn,As)As mixed crystals 8| 
with experimental data. This finding opens a question whether and why the number of the 
As antisites is correlated with the level of Mn doping. 

We investigate the correlation between the donors and acceptors in partly covalent III-V 
semiconductors such as GaAs. The cohesion energy of the covalent networks has a max- 
imum if the Fermi energy lies within a band gap (loffe-Regel's rule j^) . Whenever the 
Fermi energy is situated in the valence or conduction band the strength of the bonds is 
reduced because the unfilled bonding states or occupied antibonding states appear, respec- 
tively. It is natural to expect that this mechanism, connected with the changes of the Fermi 
energy position in dependence on the compensation, leads to some kind of acceptor - donor 
correlation. 



II. CORRELATION ENERGY 



We consider a crystal doped with both acceptors and donors. For simplicity, we assume 
that both acceptors and donors are substitutional impurities, which is the case of both As 
antisites and Mn in GaAs. Their concentrations are xa and xn, respectively. The total 
energy of the doped crystal, normalized to a unit cell, is W{xa,X£)). We show first that 
the derivatives of W{xa,xd) with respect to the concentrations xa and xd determine the 
formation energies of the defects, and also the correlation energy of the co-doping. 

We start with a large unit cell (LUC) consisting of N unit cells of the mixed crystal. The 
formation energy of an acceptor A is defined as the reaction energy of the substitution 
process 

LUC(^) + A — > LUC(^) + X. 
Here, LUC^^) is a large unit cell with one extra acceptor A replacing an atom X of the 
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oriEf inal LUC(^). In our notation, the corresponding reaction energy is 

Xd)=N- {W{xa + 1/iV, Xd) - W{XA, Xd)) + Eaton.{X) - Eatom{A). (1) 

The last two terms in Eq. (1) are the total energies of free-standing atoms X and A, respec- 
tively. The additional constant Eatom{.X) — Eatomi^) , though crucial for the correct absolute 
value of the formation energy, does not depend on the actual composition of the material 
and so it is not important for the concentration-dependent effects we have in mind. With 
increasing size of the large unit cell, — > oo, the first term in Eq. (1) approaches the 
derivative of W{xa,xd) with respect to xa- We have 

E\XA, xd) = ^^^g^^''''^ + EatoUX) " EatoU^), (2) 

in a close analogy with the definition of the chemical potential. Similarly, the formation 
energy of a donor D substituting for an atom Y is 
T?Df N dW{xA,XD) 

E [Xa, Xd) = h Eatom{Y) - Eatom[D), (3) 

Finally, the compositional dependence of the formation energies can be characterized by the 
correlation energy 

K(xA xd) = ^ = ^ = (4) 
dxD dxA dxAdxD 

The correlation energy can be alternatively, for finite large unit cells, expressed in terms 
of the four total energies corresponding to the reference LUC^'^^^ and to related systems 
with an extra acceptor (LUC^'^^-'), extra donor (LUC*^^^^), and both acceptor and donor 
(LUC(^^)), 

The correlation energy K{xa, xd) is positive if the formation energy of one impurity increases 
in the presence of the other. This means, in the case of semiconductors, that the material 
tends to be either n-type or p-type rather than a compensated semiconductor. On the other 
hand, negative correlation energy indicates that the presence of impurities of one kind makes 
the incorporation of the other dopants easier. In this case we can speak about a preferential 
compensation. 

To investigate the correlation among the dopants, we use the coherent potential approx- 
imation (CPA) [kJ. There are two reasons for this. First of all, CPA describes the config- 
urationally averaged behavior of the mixed systems rather than particular arrangements of 
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Figure 1: Dependence of the As antisite formation energy on the concentration x of Mn acceptors 
for two series, Gai_2,.Mni.As (boxes, y = 0.00) and Gao.99_a;Mn2,.Aso.oiAs (circles, y = 0.01). The 
changes /S.E of the formation energy with respect to the reference systems with x = 0.04 are shown. 

the impurities. This fits well to the thermodynamical Yoffe-Regel mechanism. In addition, 
the CPA is applicable to impure crystals with arbitrary concentrations of the dopants and it 
is particularly suitable for the description of the compositional dependence of the electronic 
properties. 

III. RESULTS 

We considered a series of multicomponent mixed crystals Gai-aj-j^A^jDj^As with acceptors 
A=Mn, Zn and donors D=As, Sn. The highly degenerate nonmagnetic mixed crystals with 
Zn were studied in parallel to the diluted magnetic semiconductor (Ga,Mn)As to make 
clear that the acceptor-donor correlation does not depend on the presence of the magnetic 
moments. With our choice, all substitutions take place in the cationic sublattice. Only 
doping giving p-type materials was considered. 

The ab initio TB-LMTO version of the CPA [11] was used to calculate the electronic 
structure of the impure crystals and of their total energies. To obtain the derivatives of the 
total energy with respect to the concentrations, we varied both x and y on a fine mesh with 
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5x, 5y = 0.005. 

Table I summarizes the calculated correlation energies K{x,0). The correlation energy 
of the co-doping is negative in all considered materials. The absolute value of K is of order 
of a few electronvolts and it generally decreases with increasing level of the doping. The 
weakest correlation if found in (Ga,Mn)As compensated with Sn. The correlation energies 
for Sn are approximately two times smaller than the correlation energies for the As antisite, 
which is a double donor, as expected. 

TABLE I 

Correlation energy for acceptors (A) and donors (D) in p-type mixed crystals Gai_a;Aa;As 
at various levels of p-type doping. 



X 


A=Mn, D=As 


A=Zn, D=As 


A=Mn, D=Sn 


0.04 


-3.36 eV 


-5.48 eV 


-1.39 eV 


0.05 


-2.67 eV 


-4.67 eV 


-1.22 eV 


0.06 


-2.21 eV 


-4.19 eV 


-1.09 eV 



The negative correlation energy of the co-doping means that the formation energies of 
both acceptors and donors decrease in the presence of the compensating impurities. This is 
shown explicitely for the substitutional Mn and for the As antisite defect in Figs. 1 and 2. 
To avoid the technical problem with the additional constants in Eqs. (2,3), we do not plot 
the entire formation energies 2 eV), but only their changes with respect to the reference 
material, mostly Gao.96Mno.04As. 

Fig. 1 shows how the formation energy of the As antisite defect changes with the con- 
centration of Mn. In addition to the main series of the data for Gai_a;Mn-rAs we considered 
also a series of mixed crystals Gao.99_a;Mnj.Aso.oiAs already containing a small portion of 
the As antisites. In both cases, formation energy is reduced by approx. 0.1 eV if the Mn 
concentration increases by a few atomic percent. 

The formation energy of the substitutional Mn was calculated for two series of mixed 
crystals, Gao.96-2/Mno.o4Asj,As and Gao.94_2^Mno.o6Asj/As. Its variation with increasing am- 
mount of the As antisites is shown in Fig. 2. The changes of the formation energy are again 
of order of 0.1 eV. 

To explain the results we turn to a simple model. In the one-electron picture, the total 
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Figure 2: Change of formation energy, AE, of substitutional Mn in Gao.96-j,Mno.o4AsjyAs (boxes) 
and Gao.94-yMno.o6-A-SyAs (circles) due to the As antisite defects as a function of their concentra- 
tion y. 

energy can be expressed in terms of the density of states g{E), 



The dependence of W{xa,xd) on the chemical composition arises both from the redistribu- 
tion of the electron states in the valence band due to the impurities and from the changes of 
the position of the Fermi level. In the case that giEp) does not change much with the chemi- 
cal composition, the total energy depends on the concentrations of the impurities mostly via 
the position of the Fermi level. Because the extra acceptors (donors) push the Fermi level to 
lower (higher) energies, the correlation energy for the acceptor-donor co-doping is expected, 
according to Eqs. (5) and (6), to be always negative. Morever, as the variations of the Ep 
are inversely proportional to g{EF), one can expect that also W ■ g{Ep) ^ const.. This is 
in a good agreement with the calculated decrease of the correlation energy with increasing 
concentration of Mn (cf . Table I) . 




(6) 
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IV. SUMMARY 

We showed that the formation energy of As antisite defects in GaAs decreases with 
increasing concentration of Mn or Zn in the cationic sublattice. The formation energy is 
reduced by approx. 0.1 eV so that the number of these native defects can be largely enhanced 
in the presence of Mn. This effect may contribute to the self-compensation behavior of 
(Ga,Mn)As mixed crystals. 

At the same time, the formation energy of the substitutional Mn is similarly reduced in 
the presence of the As antisites or by intentional co-doping with Sn. This can be interpreted 
that the presence of the donors, either native defects or intentional dopants, is important 
for an improved solubility of Mn in III-V materials. 

Finally, the correlation between donors and acceptors seems to be a general feature 
tending to a selfcompensation. 
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